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Using a numerical library for arbitrary precision arithmetic I study the 
irregular dependence of the diffusion coefficient on the slope of a piecewise 
linear map defining a dynamical system. I find that the graph of the diffu- 
sion coefficient as a function of the slope has the fractal dimension I, but 
the convergence to this limit is slowed down by logarithmic corrections. 
The exponent controlling this correction depends on the slope and is either 
1 or 2 depending on existence and properties of a Markov partition. 

PACS numbers: 05.45.Pq, 05.45.Ac, 45.30,+s 



1. Introduction 

Determining the transport coefficients of many-particle systems is one of 
the fundamental problems of nonequilibrium statistical physics. This is also 
a notoriously difficult problem: it turns out that even in a simplified system 
where a single particle moves in a periodic array of scatterers, the drift and 
diffusion coefficients are highly irregular, apparently nowhere differentiable 
functions of control parameters ^ [21 El H] (an example is also shown in 
Fig- [U below). Closer inspection of these functions reveals that usually their 
"irregularities" are not random, but rather form patterns. This has led some 
researchers to the idea that the graphs of these functions are fractals with 
nontrivial, perhaps locally varying fractal dimension. This is an interesting 
concept, since fractal structures are quite commonly found in dynamical 
systems in various contexts [21 0] • 

Until very recently no reliable investigation of fractal properties of trans- 
port coefficients was possible because all general methods of calculating 
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transport coefficients in dynamical systems - e.g. the transition matrix tech- 
nique combined with the escape rate formalism ^ |^1 EH EJ El E] , the Green- 
Kubo formula ^EllI], or the periodic-orbit formalism [31 [7j - eventually 
lead to complicated and time-consuming numerical calculations. Moreover, 
usually they are applicable only for some special values of the control param- 
eters. For these reasons none of them could be used to collect a sufficiently 
large number (counted at least in millions) of very precise data required 
in fractal analysis. This situation changed when Groeneveld and Klages |H] 
gave exact formulas for the transport coefficients in a simple one-dimensional 
dynamical system introduced by Grossmann and Fujisaka [S]. 

The system investigated by Groeneveld and Klages can be considered as 
a model of a particle moving in a one-dimensional array of scatterers. The 
role of the equation of motion is played by a one-dimensional map M 

X n+1 = M atb (x n ) (1) 

where n is a discrete-time variable and the map M a b : R — > R is given by a 
simple linear function 

M a>b (x) = ax + b, x£l (2) 

on the fundamental interval Iq = [— |, |), with a > 1, b € R being the con- 
trol parameters representing the slope and the bias of the map, respectively; 
the map is then continued periodically onto the real line by a lift of degree 
one, i.e. by requiring that 

M a , b (x + 1) = M a>b (x) + 1, x € R. (3) 

For a > 1 the Lyapunov exponent of this system is positive, and so the dy- 
namics defined by M ab is chaotic. If we choose an arbitrary number xo as 
the starting point, the resulting sequence (x n ) will almost always look "ran- 
dom" (the set of points xq generating a regular, periodic or quasi-periodic 
sequence (x n ) is of Lebesgue measure 0). For suitably chosen values of a 
and b the deterministic dynamics defined by the map M a b is equivalent to 
a Markov stochastic process of random-walk type: each deterministic tra- 
jectory {x n ) is equivalent to some particular realization of the correspond- 
ing random-walk process, and taking the average over all initial states is 
equivalent to calculating averages over the corresponding Gibbs ensamble. 
Actually any random-walk Markov process in a one-dimensional periodic 
system with fixed transition rates can be translated into the language of 
simple piecewise linear deterministic maps ^Uj. From this point of view 
there is no surprise that the process defined by M a ^ b (or similar maps) is 
called "deterministic diffusion" and that the two basic transport coefficients, 
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Fig. 1. The diffusion coefficient D for the map M a ^ as a function of the slope a 
for the bias 6 = 0. The inset depicts a blow-up of a region 6 < a < 6.0001 (pointed 
at by the arrow). 



the drift velocity J and diffusion constant D, can be defined as 

j= lim <__> D = lim ^ ~ (Xn)2 (4) 

n^oo 77, n.^oo 2fl 

where (• • • ) denotes the average over the uniform ensamble of initial values 
x . 

The graph of the diffusion coefficient D for the map M a ^ as a function 
of the slope a for the bias b = is shown in Fig. ^ The fractal properties 
of this highly irregular graph (as well as that of the drift velocity J) were 
recently studied by Klages and Klaufi Using two numerical methods: 
the box counting and the autocorrelation function methods, they found 
that the local fractal dimensions of graphs of J and D are well-defined, but 
highly irregular functions of a and b. In other words they found that the 
graph shown in Fig. ^ cannot be described with a single fractal dimension, 
but rather by a set of quickly varying local fractal dimensions. Taking 
this into account they suggested that the local fractal dimensions of graphs 
of J and D as functions of the slope a with the bias b fixed are fractal 
themselves. This, in turn, leads to the concept of a "fractal fractal dimension 
of deterministic transport coefficients" jllj . 

Klages and Klaufi found the local fractal dimension T> of J(a, b) and 
D(a,b) to be very close to 1. The autocorrelation function method gave 
1.05 < T> < 1.15, and the box-counting method gave even lower values 
1.02 < T> < 1.05. All these values are very close to 1, i.e. the fractal 
dimension of a regular, "smooth" curve. This suggests that perhaps T> = 1 
and that the results obtained in the above-mentioned computer simulations 
reflected extremely slow convergence of T> to its true asymptotic limit of 
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infinitesimally small "boxes". Such slow convergence is often caused by 
logarithmic corrections. The main purpose of my paper is thus to try and 
find such tiny corrections. To accomplish this task I will use a different 
method of calculating the local fractal dimension of a curve - the so called 
oscillation method ^2] ~~ an d I will carry out all calculations with a help of 
a special high-accuracy numerical library, which will enable me to consider 
exceptionally small "boxes" . For sake of simplicity I shall restrict the present 
study to the simplest case of zero bias (b = 0) where, due to symmetry, the 
drift velocity J = 0. 

The structure of the paper is as follows. Section |2] briefly describes 
the method I have used to determine the local fractal dimension. It de- 
scribes all steps necessary to calculate the diffusion coefficient for the map 
M a fi, a method of calculating the local fractal dimension of a continuous 
curve ("the oscillation method"), mathematical formulation of the main 
conjecture about the logarithmic convergence of the fractal dimension to its 
limiting value, and the numerical aspects of the algorithms used. Section 
OH presents the main results. Finally, section 0] is devoted to discussion of 
results. 



2. Method 

2.1. Transport coefficient 

The explicit formulas 8 for the transport coefficients in the model de- 
pend on some auxiliary variables. For each a > 1,6 € R, and e = ± we 
define two infinite sequences, (y^), r = 0, 1, . . ., and (n^.), r = 1, 2, . . ., con- 
sisting of real and integer numbers, respectively. Their values are uniquely 
determined by demanding that t/q = | and that for each r > 

n € r + y € r = ay € r _ x + b, (5) 



with additional conditions n% € 7L and y% 6 Iq, where 1$ = (— 5, 5] and 

re "N- numbers": 

r 

K = - e - + Y.<i ( 6 ) 



I = [— 2, 2)- Next we define "N- numbers" 



1 00 

r=0 

where r, k, I > 0. The basic transport coefficients, J and D, can be now 
expressed (jHj, cf ^Hj) as 

, _ N 2 ,o jV 3 ,o ~N 2:1 J + N h2 J 2 
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2.2. The oscillation method 

The fractal dimension T> of a continuous function / defined on an interval 
[io,£i] can be evaluated through analysis of its Holder exponents jT2j and 
r-oscillations ^2j. The r-oscillation of / at t G (to,ti) is defined as 

osc T (/;t)= sup /(*')- mf /(*')■ (9) 

\t-t'\<T \t-t'\<T 

If there exist constants c > and < H < 1 such that for all r 

osc T (t) < ct h (10) 

then / is called a Holderian of exponent H at t and its fractal dimension P 
at i is related to H through 

V<2-H. (11) 
Similarly, if there exist constants c > and < H < 1 such that for all r 

osc T (t) > ct h (12) 
then / is called an anti-Holderian of exponent H at t and 

V>2-H. (13) 

£.5. Reformulation of the problem and numerical implementation 

We are now ready to formulate our main conjecture: for the map M a j, 
with 6 = 

° SCT(L>;a) w c(a) [- log(r)]^ a ) , as r - 0, (14) 
r 

where the prefactor c(a) > and the exponent 7(a) > 0. Owing to © - 
(|13|) this implies that the fractal dimension of the graph of D as a function of 
the control parameter a is equal 1, but the convergence to this limiting value 
is logarithmically slow, with the exponent j(t) controlling the convergence 
rate. 

I checked conjecture (|14|) using equations © - ©• Because the consecu- 
tive terms of the sequences nf in © can be generated extremely efficiently 
(it took less than a second to generate 20 000 data points used to draw 
Figure I decided to employ GMP (GNU multiple precision arithmetic 
library, version 4.1.2) - a library for arbitrary precision arithmetic [Tl]. Al- 
though any calculations performed by such a library must be several orders 
of magnitude slower that those performed directly, thanks to the GMP I was 
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able to study numerically the limit of r —* for r ranging from 1 down to at 
least 10 -100 (the accuracy of typical computers currently available, without 
using such special-purpose libraries, is limited typically to about 10~ 17 ). 
Such fine resolution of my calculations will turn out crucial for determining 
logarithmic corrections predicted by formula l|14|l. 

3. Results 

In my present study I decided to concentrate on verifying conjecture 
(j!4j) for the symmetrical case of vanishing bias (b = 0) and for several se- 
lected values of the slope a. I decided to examine especially carefully the 
slopes corresponding to finite Markov partitions. These values are partic- 
ularly interesting from the theoretical point of view for example, 
even though the diffusion coefficient D is always continuous in a [S], the 
numerator and denominator in Q can be discontinuous for such (and only 
such) slopes. One of the reasons why I chose to set b = is that in this 
case the slopes corresponding to Markov partitions can be found easily on 
a computer: they are algebraic numbers generating periodic sequences rik 
and yk. Such values of a will be henceforth referred to as "Markov slopes". 

Figure [21 presents results obtained for two Markov slopes generating 
strictly periodic sequences (i.e., y e k+l = y\ for a period I > and all k > 0). 
The value used in Fig. [2^ is a = 3. This slope generates simple sequences y k 
and yT of period 1. The oscillations of the diffusion coefficient on intervals 
(3 — |r, 3 + \t), rescaled by r, obtained for r = 10~ n ,n = 1 . . . 100, are 
represented in this plot by circles, while the solid line represents a quadratic 
fit of form px 2 + qx + r with p = 0.733, q = 0.57, and r = —0.2. The fit is 
excellent, which suggests that for this slope the exponent 7, as defined in 
eq. (|14|). is 2. I have obtained similarly good quadratic fits for other integer 
odd slopes, which also generate simple sequences of period 1 0. 

A question arises whether the length of the period has any influence on 
the limiting value of osc r (P; a)/r. FigureEb shows the results obtained for 
the largest root of the polynomial a 4 — 4a 3 + 2a 2 + 3, i.e. for a ~ 3.314. For 
this slope the period I = 4. Just as in the previous example, osc r (P;a)/r 
can be approximated by a quadratic, although the fit is not as excellent as 
for a = 3. The inset in this figure depicts the blow-up of the data obtained 
for 1 > r > 10~ 20 . It shows that for 1 > r > 10~ 10 the value of osc r (£>; a)/r 
fluctuates about a constant value, which might suggest that the exponent 7 
vanishes. This example demonstrates that the "resolution" of calculations 
used in ref. e.i. r m ; n < 10~ 10 , is insufficient to find the true asymptotic 
fractal properties of the graph of the diffusion coefficient as a function of 
the slope a. 

All other Markov slopes generate sequences n| and yi with the period 
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starting at a term m > 0. There are two cases: either all terms of y£ 
are different from all the terms of y^ or the two sequences have the same 
period sequence. Let's start from the latter case. Three examples of results 
obtained for Markov slopes of this type are shown in Figure The data 
depicted in this graph were calculated for a = 4 (circles), for the largest 
root of a 4 — 4a 3 + 2a 2 + 2, i.e. for a « 3.3497 (squares), and for the largest 
root of a 5 - 4a 4 + 2a 3 + 3a 2 - 2a + 4, i.e. for a 3.043 (pluses) . For the first 
two of the above slopes all the terms of the sequences yt eventually vanish. 
In the case of a = 4 this happens for k > 1, and for a ~ 3.3497 the terms of 
y e k vanish for k > 4. As for the third example, a ~ 3.043, the period starts 
at k = 2 and has the length 6. Specifically, the first five terms of y~£ are 
(to 3 S.F.) 1/2, -0.478, -0.456, -0.387, -0.179, and for k > 5 the terms 
satisfy yt = —yi_% = 2//T_ 3 - As we can see, in all these cases osc T (2?; a)/r 
diverges linearly with log(r _1 ), suggesting that conjecture (|14|) is satisfied 
with the exponent 7 = 1. 

The results obtained for the last category of Markov slopes, i.e. slopes 
generating two disjoint periodic sequences y^ and y^ , are shown in Fig. 03. 
The two data sets were collected for the largest root of a 5 — 4a 4 + 2a 3 + a 2 + 
4a— 2, i.e. for a ~ 3.16 (circles) and for the largest root of a 5 — 4a 4 — 4a 2 +a+4, 
i.e. for a « 4.20 (squares). The first of these slopes, a 3.16, generates a 
sequence with a period of length 3 starting at y^, while a ~ 4.20 generates a 
sequence of period length 4 starting at yf. As we can see, for these Markov 
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Fig. 2. Oscillations of the diffusion coefficient rescaled in accordance with eq. I|14|) 
calculated on intervals (a— \t, a+^r) for r = 10~™, n = 1, . . . , 100 and for (a) a = 3 
and (b) a « 3.314. The circles represent results calculated from conjecture (|14J) . 
and the solid lines are quadratic fits. The inset in graph (b) shows the blow-up of 
results obtained for 1 > r > 10~ 20 . 
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slopes osc r (P; a)/r can be very well approximated by a quadratic function 
of log(r _1 ). This suggests that in this case conjecture (fTl)) is satisfied with 
the exponent 7 = 2. 

Figure 0] presents results for two slopes a which do not correspond to a 
Markov partition. These are a = it ~ 3.14 (left panel) and a = 9/2 (right 
panel). The results turned out to be "noisy" and I had to increase the 
resolution down to r m i n = 1CT 300 . As can be seen, in both cases osc r (P; o)/t 
has a clear trend linear in log(r _1 ). However, this limiting behaviour is 
disturbed by very large fluctuations. Actually, for a = tt these fluctuations 
are so large that my numerical results cannot rule out the possibility that 
in this case the limit lim r ^o osc r (P; a)/r does not exist. 



4. Conclusions 

In my study I proposed and verified numerically a hypothesis that the lo- 
cal fractal dimension of the graph of the diffusion coefficient of the piecewise 
linear map M a ^ as a function of the slope a is T> = 1 and that the conver- 
gence to this limit is slowed down by logarithmic corrections described by 
eq. (|14[) . This contradicts the earlier findings of Klages and Klaufl that 
this graph is a fractal with a locally varying fractal dimension T>(a) > 1 
which, when plotted as a function of the slope, forms a fractal itself. 

I found that the exponent 7, which controls the logarithmic correction, 
is actually a function of the slope a. Interestingly, 7 appears to be a dis- 
continuous function that can take only one of two values: 1 or 2. The value 




-log(T) -log(x) 



Fig. 3. The same as in Fig. El (a) a = 4 (circles), a w 3.3497 (squares), a « 3.043 
(pluses) with linear fits (solid lines); (b) a « 3.160 (circles), a w 4.200 (squares) 
with quadratic fits (solid lines). See text for details. 
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of 1 corresponds to Markov slopes that generate two disjoint sequences y^ 
and y^ , and the value of 7 = 2 corresponds to Markov slopes that gener- 
ate periodic sequences y^ and y^ with the same period terms. The case 
of slopes that do not generate Markov partitions is not clear - apparently 
7 = 1, but this statement cannot be verified numerically because of very 
large fluctuations making the convergence extremely slow. Note that my 
findings imply that both of the sets {a: 7(a) = 1} and {a: 7(a) = 2} are 
dense, and hence that 7(a) is nowhere continuous. 

Any numerical study is naturally restricted to investigation of several 
particular cases. It cannot be ruled out that I have overlooked some cate- 
gories of slopes a for which the local fractal dimension of the graph of the 
diffusion coefficient behaves in a way not predicted in this study. Numerical 
investigation of the logarithmic correction to the asymptotic limit is here a 
particularly delicate problem, as apparently 7 is a nowhere continuous func- 
tion of the slope. Only through an analytical approach could this problem 
be solved conclusively. Such a study will be published elsewhere. 
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